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We present the results of a radial velocity (RV) survey of 14 brown dwarfs (BDs) and 
very low-mass (VLM) stars in the Upper Scorpius OB association (UScoOB) and 3 BD can- 
didates in the p Ophiuchi dark cloud core. We obtained high-resolution echelle spectra at the 
Very Large Telescope using Ultraviolet and Visual Echelle Spectrograph (UVES) at two dif- 
ferent epochs for each object, and measured the shifts in their RVs to identify candidates for 
binary/multiple systems in the sample. The average time separation of the RV measurements 
is 21.6 d, and our survey is sensitive to the binaries with separation < 0.1 au. We found 
that 4 out of 17 objects (or 24l]J per cent by fraction) show a significant RV change in 4— 
33 d time scale, and are considered as binary/multiple 'candidates.' We found no double-lined 
spectroscopic binaries in our sample, based on the shape of cross-correlation curves. The RV 
dispersion of the objects in UScoOB is found to be very similar to that of the BD and VLM 
stars in Chamaeleon I (Cha I). We also found the distribution of the mean rotational velocities 
(v smi) of the UScoOB objects is similar to that of the Cha I, but the dispersion of v sini is 
much larger than that of the Cha I objects. 



Key words: stars: binaries:spectroscopic - 
- stars:planetary system: formation 



stars: low-mass, brown dwarfs - stars: formation 



1 INTRODUCTION 

Most stars are member of binary systems and it is therefore im- 
portant that a complete star formation theory be able to predict 
the binary fraction, period distribution, and mass-ratio distribution 
of stellar objects across a wide range of masses. Furthermore, the 
study of individual binary systems is the only direct means to de- 
termine fundamental stellar properties such as stellar masses and 
radii. 

Recent high-resolution imaging studies of brown dwarfs 
(BDs) and very low-mass (VLM) stars have placed strong con- 
straints on binaries with separations of ~ 1 — 100 au. For exam- 
ple, Hubble Space Telescope (HST) observations of a Per and the 
Pleiades indicates a binary fraction (/) of > 10 per cent with a bias 
towards separations (a) of less than 15 au, and a mass-ratio (g) of > 
0.7 iMartfn et all2003 | ) for objects arou nd and below the hydrogen 
burning limit (see also lBouv et alj2006t) . A similar lack of wide bi - 
naries was found in the field T-dwarf study iBurgasser et all2003h . 
while / « 15 per cent (for the fiel d objects with 0.03 -0.1 Mq 
and M8.0-L0.5) was determined bv lClose etafl <2003l) using the 
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adaptive optics at Gemini North. They also found the vast majority 
of sy stems have a semimajor a xis < 20 au (se e also lSiegler et alj 
120051 but note too the results of lLuhmanll2004. An HST study of 
more than 80 field late M and L dwarfs iGizis et alJ2003l) indicated 
/ ~ 15 per cent with separations in the range of 1.6-16 au. For a 
small (12) sample of BDs and VLM stars (0.04-0.1 M fl ) in Upper 
Scorpi us OB association (UScoOB'l. lKraus. White. & Hillenbrand! 
( 2005) found / = 25lg 6 per cent for 5 au < a < 18 au by using a 
similar imaging technique. More recently. iBasri & Reinersl feoOrJ) . 
combined with the results of earlier works, found the upper limit of 
the overall binary fraction for VLM stars of 26 ± 10 per cent. 



Using a Monte Carlo simulation, the data from radial veloc- 
ity surveys available in the literature, and by carefully consider- 
ing the sensitivity and sampling biases, iMaxted & Jeffries! 120051) 
found an overall BD/VLM binary frequency of 32^15 per cent as- 
suming f = 15 per c ent for a > 2.6 au. A recent photometric study 
IPinfield et all200 3) of low-mass objects in Pleiades and Praesepe 
suggested, albeit indirectly, / as large as 50 per cent, which would 
only be compatible with direct imaging studies if 70-80 per cent 
of those binaries have a < 1 au. For a more comprehensive review 
of the current status of BD/VLM binary fraction and the separa- 
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tion distri bution, readers are refer to a recent review of multiplicity 
studies bv lBurgasser et allfcOOd) . 

The extensive imaging surveys provide excellent observa- 
tional constraints on wider BD+BD binaries, but it is now nec- 
essary to search for shorter period BD+BD binaries systemati- 
cally. Binaries with the separation of less than 1 au are not re- 
solved by current imaging techniques, but will be detectable as 
spectroscopic binaries, providing the mass-ratio is not too extreme, 
and velocity separation is large enough. The first BD+BD spec- 
troscopic binary, PP1 15 teasri & Maitfill999h . showed a double- 
peaked cross-correlation function with a maximum velocity sep- 
aration of > 70kms _1 . The binary was found to have an ec- 
centric orbit (e = 0.4) with a period of ~ 5.8 d. lBasri & Martini 
1 19991) suggested that the formation process of substellar objects 
is biased towards smaller separation binaries based on the short 
period of PP1 15 and the lack of Pleiades BD binaries with separa- 
tions > 40 au. Note that the median separation of binarie s with 
solar-type primaries is 30 au iDuquennov & Mavodll99ll) . Pio- 
neering work on the RVs o f BDs and VLMs are pre sented by 
iGuenther & Wuchterll J2003h . iKenvon et al] fc005l) and Ijoergensl 
(2006lJ who found a several binary candidates; however, the or- 
bital parameters and masses of binaries remains unknown because 
the follow-up spectroscopic monitoring is lacking or is still being 
undertaken. In addition to the follow-up observations, the number 
of BDs and VLM star binary candidates needs to be increased in 
order to have better statistics on short-period binary parameters. 

The first BD+BD ecli psing binary (2MASS J0532184- 
J0546085) was discovered bv lStassun. Mathieu. & Valentil feood) 
from the J-band photometric monitoring of the system. Combining 
their light curves and the results from the follow-up radial velocity 
measurements, they were able to determined the precise orbital and 
physical parameters of the system. The projected semimajor axis 
and the period of the binary are found as 0.0398 ± 0.0010 au and 
9.779621 ± 0.000042 d respectively. 

The separation distribution of BD/VLM binaries is criti- 
cal to understanding their origin. There are two main mod- 
els for the formation of BDs and VLM stars: first, they have 
low masses be cause they form in low-m ass, dense molecular 
cloud cores (e.g. lPadoan & Nordlunj2002l) : second, BD/VLM ob- 
jects have low masses because they are ejected from the dense 
core in which they form via dynamical interactions in multi- 
ple system, cu tting off their accretion before they have reached 
stellar masses iReipurth & Clarkdl200ll : lBate. Bonnell. & Bromra 
l2002ah . Alternatively, there is a third model in which a free- 
floating BD or planetary-mass object can be formed in the pro- 
cess of the photo-evaporation (e. g. iMcCaughrean & Andersenl 
l2002tlWhitworfh & Zinneckej2004h with the outer layers of a pre- 
stellar core (~ 0.2 Mq) removed by the strong radiation pressure 
from the nearby massive OB stars before the accretion onto the 
protostar at core centre occurs. 

Due to the dynamical interaction involved in the second 
model, BD/VLM binaries that survive are generally expected to 
have small separations. In the first model, wider binaries may be 
expected to be more common. iBate. Bonnell. & Bromml l2002bl) 
suggested that close binaries (a < 10 au) do not form directly, 
but result from hardening of wider systems though a combination 
of dynamical interactions, accretion and interactions with circum- 
binary discs. If BD/VLM binaries have formed through such mech- 
anisms, one would not expect to find binaries with 1-10 au sepa- 
rations without also finding many with separation < 1 au. If an 
absence/rarity of binaries with 1 au were found, it may support the 



idea that they are ejected quickly from multiple systems before they 
have undergone the interactions that shorten their periods. 

Our immediate aim is to identify spectroscopic and close 
BD/VLM binaries using the high-resolution echelle spectroscopy 
at two epochs. This experiment is sensitive to VLM binaries with 
separations of < 0.1 au which corresponds to a period of ~ 10 d. 
A larger sample of candidates will enable us to measure the binary 
fraction of these short-period/close binaries (once confirmed), and 
address whether there is a significant population of 'hidden' VLM 
companions. The long term goal of this project to follow up the 
binary candidates found in this paper by spectroscopically mon- 
itoring them over different time scales, enabling us to obtain the 
radial velocity curves and their minimum masses. 

In Section|2| we describe the observations and the data reduc- 
tion. The results of radial velocity and rotational velocity (v sin i) 
measurements are presented in Section [3] We discuss the bi- 
nary/multiplicity fraction indicated by our RV survey in Section|4] 
and give our conclusions in Section[5] 



2 OBSERVATIONS 

Our sample consists of 18 young, very low-mass o bjects: 15 objects 
in the Upper Scorpius OB association (d ~ 145 pc.lde Zeeuw et alj 
1 19991) from the list of lArdila. Martin. & Basril <2000l) and 3 ob- 
jects in the p Ophiuchi cloud core (d w 150 pc, Ide Zeeuw et alj 
1 19971) from lLuhman & Riekd Jl999h . The spectral type of the ob- 
jects range between M5 and M8.5, and the age <~ lOMyr 
jLuhman & Riekdl 1 9991 : lArdila et all200(i iMuzerolle et all2003t 
Kraus et al. 2005). The sample is not complete, and the selection 
was solely based on brightness and the observability. The basic 
properties of the targets based on the literature is summarised in 
TableQ 

We obtained high-resolution spectra with the Kueyen tele- 
scope of the Very Large Telescope (VLT-Cerro Parnal, Chile) using 
the UVES echelle spectrograph. The observations were carried out 
between 2004 April 5 and 2004 May 17 in service mode. For each 
object, spectra were obtained at two different epochs separated by 
4-33 d. For each object on a given night, two separate spectra were 
obtained consecutively. This allows us to derive more re liable un- 
certain ty estimates in the RV values of our targets (c.f. Ijoergensl 
l2006bh ." The data were obtained using the red arm of UVES spec- 
trograph with two mosaic CCDs (EEV + MIT/LL with 2kx4k pix- 
els). The wavelength coverage of 6708 - 10,250 A and the spectral 
resolution R ~ 40, 000 were used. The slit width and length of 1" 
and 12" were used respectively with a typical seeing of 0.8". 

The data were reduced via the standard ESO pipeline proce- 
dures for UVES echelle spectra. In summary, the data were cor- 
rected for bias, interorder background, sky background, sky emis- 
sion lines and cosmic ray hits. They were then flattened, optimally 
extracted, and finally the different orders were merged. No bin- 
ning was performed to achieve high resolution required for the RV 
measurements. The wavelength was calibrated using the Thorium- 
Argon arc spectra with a typical value of the standard deviation of 
the dispersion solution of 5 mA which corresponds to 0.2 km s _1 
at the central wavelength 8600 A. However, the autoguiding of the 
telescope keeps the star at the centre of the slit with about a tenth of 
the FWHM (1 km s _1 ) which sets the upper limit for the system- 
atic error in the RV measurements IB ailer- Jone J2004T) . In fact, we 
find that any systematic error is negligible compared to the random 
error associated with our radial velocity measurements (see Sec- 
tion !3. II . A typical signal-to-noise ratio (S/N) per wavelength bin 
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Table 1. Summary of known properties of the targets from litera tur e: a.ETuhman & Riekd <1999l) (original l i st for p Ovh). b.lArdila. Martin. & BasrilfcOOC 
^origin al list for UScoOB), c. lWilking. Greene. & MeveJ ll999t) . d. iMuzerolle et alj <2003l) . e. lKraus et all l2005t) . and f. iMohantv. Javawardhana. & Basi 



Object Sp. mass RV v sin i Known multiple? 

[M Q ] [kms" 1 ] [kms" 1 ] 



GY 5 


M7 C 


0.07 d 


-6.3 ± 1.9 d 


16.8 ± 2. 7 d 


no 


GY 141 


M8.5 a 


0.02 d 




6.0-f 


no 


GY310 


M8.5 C 


0.08 a > d 




10.0-f 


no 


USco 40 


M5 b 


0.1 b 




37.5* 


no 


USco 53 


M5 b 


0.1 b 




45.0-f 


no 


USco 55 


M5.5 b 


0.10 + 0.07 c 




12.0-f 


visual 


USco 66 


M6 b 


0.07 + 0.07 c 


-4.4±0.6 d 


27.5-f 


visual 


USco 67 


M5.5 b 


0.10 c 




18.0 f 


no 


USco 75 


M6 b 


0.07 c 


-5.6 ± l.l d 


63.0 f 


no 


USco 100 


M7 b 


0.05 e 


-8.9±0.6 d 


50.0-f 


no 


USco 101 


M5 b 


0.05 b 






no 


USco 104 


M5 b 


0.05 b 




16.0^ 


no 


USco 109 


M6 b 


0.07 + 0.04 c 


-3.8±0.7 d 


6.0^ 


visual 


USco 112 


M5.5 b 


0.1 c 




8.0 f 


no 


USco 121 


M6 b 


0.02 b 


-38.9 ± 1.0 d 




no 


USco 128 


M7 b 


0.05 e 


-3.0 ± 1.6 d 


0.0^ 


no 


USco 130 


M7.5 e 


0.04 e 




U.0 f 


no 


USco 132 


M7 b 


0.05 e 


-8.2 ± l.l d 




no 



of the spectra is about 15, and the heliocentric velocity correction 
was applied to the final spectra. 

In addition to the main targets in TableQ we also obtained the 
spectra of LHS 49 (on 2004-April-17) and HD 140538 (on 2004- 
May-06), which are the RV template and the RV standard stars re- 
spectively. These data were obtained using the same telescope and 
the same instrument setup as for our main targets. The S/N of the 
both objects were about 60. 



3 RESULTS 

3.1 Radial velocities 

The radial velocities of each object were determined by using 
the cross-correlation function of the object spectrum with that of 
a template star which has a similar spectra type. By visual in- 
spection, the wavelength ranges used for the cross-correlation cal- 
culations are chosen by avoiding the regions of spectra affected 
by the telluric lines, and defects and fringes (in near infrared) 
of the CCDs. The radial velocities of objects with respect to the 
template are obtained by measuring the location of the peak in 
the cross-correlation function. The location of the peak is deter- 
mined by fitting the section of the cross-correlation function around 
the peak by a second order polynomial. LHS 49 (Proxima Cen, 
M5.5) was chosen as the template for this purpose. The helio- 
centric radial velocity of the template object LHS 49 was ob- 
tained by measuring the wavelength shifts of the prominent photo- 
spheric absorption features Kl AA7664.911, 7698.974. This gives 
us RVlhs49 = — 22.6 ± 0.5 k m s -1 , which is in good agree ment 
with the earlier measurement of lGarcia-Sanchez et alj 1200 ll) who 
found RVlhs49 = - 21.7 ± 1.8 km s" 1 . The heliocentric RV of 
each object can be then calculated by adding RVlhs49 with the 
RV of each object with respect to LHS 49. In the following mea- 
surements of the heliocentric radial velocities, our measurement 
(RVlhs49 = — 22.6 ± 0.5 km s _1 ) will be used for consistency. 



Before applying the cross-correlation technique to our main 
targets, we have applied the technique to the radial velocity stan- 
dard HD 140538 (G2.5V) for which an high-accuracy RV mea- 
surement via th e fixed-configuration, c ross-dispersed ecelle spec- 
trograph Elodie iBaranne & et alll99ffl is available. This was done 
so to ensure not only the validity of the cross-correlation technique, 
but also the validity of the wavelength calibration. In this test, we 
found the heliocentric RVhd 140538 = 19.17±0.40kms~ 1 which 
is in good agreement with the Elodie radial velocity me asurement 
of 19.00 ±0.05 kms" 1 lUdrv. Mayor. & OueloJ 19991) . Note that 
the uncertainty in RVhdi40538 was computed by combining (in 
quadrature) the uncertainty (0.4 km s" 1 ) in the heliocentric RV of 
the template star HD 74497 (G3V) and the uncertainty in the rel- 
ative RV of HD 140538 (0.02 km s" 1 ), measured by cross corre- 
lation with the template star. Clearly, the former dominates in the 
uncertainty of RVhd 140538 shown above. The uncertainty in the 
heliocentric RV of the template is rather large because it was deter- 
mined from a few photospheric lines. 

The results of the heliocentric RV measurements (from two 
epochs for each object) are summarised in Table 2 along with their 
uncertainties. The uncertainties listed in columns 4 and 6 include 
the uncertainty of the absolute RV of the template star, added in 
quadrature with the uncertainty found from the cross correlation 
analysis. On the other hand, the uncertainties in the relative radial 
velocities (orrv) listed in column 5 of the table do not include the 
uncertainty of the template star. It is ctrrv that is relevant for the 
identification of RV variable objects. 

To determine the ctrrv, we have used the following proce- 
dure: (1) add random gaussian deviate noise to the object spectra 
using the corresponding variance spectra 1 , (2) compute the cross- 
correlation curve using the spectra with added noise and the tem- 
plate, (3) measure the RV by locating the peak, (4) repeat 1-3 for 



1 Defined by the va riance of the \ 2 fit to the signal obtained during opti- 
mum extraction (c.f. iHornel 1986ft . 
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Table 2. Summary of the observations, the heliocentric radial velocities (RV) from two-epoch and the average rotational velocities (v sini). The uncertainties 
of relative radial velocities (cr RRV ) with respect to the template star LHS 049 and the average radial velocities (RV) are also given. N.B. the uncertainties 
listed along with the heliocentric radial velocities in column 4 and 6 include the uncertainty in the heliocentric RV of the template star, but ctr rv in column 5 
do not. 



Object Date HJD-2453100 RV o- RRV RV vsini 

[kms -1 ] [kms -1 ] [kms -1 ] [kms -1 ] 



GY 5 


onn/i a o/i 
2UU4-Apr-24 


on "7 1 C\Q A A n 

2U. 1 iyo44u 


a r 7'7 _L o no 
— 0. t ( ± 2.U2 


1 Ci£. 

i.yo 












onn/i \/inir f\n 
2UU4-May-U / 


1 1 T70 i ni 


a no _i_ not; 
— 0.02 ± U.oO 


0.68 


p. Qn _i_ o i n 

— o.oy ± 2.ly 


16. 


5 ± 


.6 


( "■ v 1/11 


zuu4-May- 1 u 


ifi. A~71 o 1 no 
JO.O 1 121U8 


< if 1 n £?q 
— 4.10 ± U.OO 


n io 

KJ.DO 












onn/i iv 4 ,-, . , i T 
zUU4-JVLay- 1 / 


43.0O83UU / 


o i net; 
— O.40 ± U.OO 


n ai 
U.4J 


— o.ol ± U.yl 


4 


1 _i_ 1 
4 ± 1 


.4 


c~v i 1 n 
Lji 31U 


Onn/1 Anr O/I 

2UU4-Apr-24 


on Qni^n/i 
2U.8J /J3U4 


—o.z/ ± u.yz 


n hi 
U. / / 












2UU4-May-uy 


tc ooo/in^o 
jO.8224U02 


o i n _L n *71 

— o.iy ± u. 7 1 


U.Dl 


— 0. /O ± l.lO 


11 


1 ± 6 


() 


T TC^ii /in 
UoCO 4U 


zUU4- Apr-ID 


1. /OOl /j4 


7 Q/l J_ n c;o 
— < .o4 ± U.02 


n i ^ 












onn/i AyioT, m 
zuu4-May-u / 


j J. /02UJ 04 


fl OQ _L n c;o 
— D.Zo ± U.OZ 


n i ^ 
U. 1 


— O.ol ± U. / 4 


Q 1 
■11 


Z ± U 


r 

.0 


UoCO jj 


onn/i Anr n/i 
2UU4-Apr-U4 




— 0.2 / ± U. / U 


n ^n 












onn/i ijnn no 
2UU4-May-U2 


OO Tin/lA*7/i 

28. /jy4o /O 


V OO _L n i^Q 

— 1 .22 ± U.OO 


0.39 


— r . / o ± u.yo 


40. 


± 


.6 


USco 55 


onn/i a 
2UU4-Apr-Uo 


1 o /i o o o m 
l.o4ZZoU/ 


c oi _L n c;n 
— 0.21 ± U.oU 


0.04 












onn/i A/inTT no 
2UU4-May-U2 


oo 01/111 no 
28.8141 iyo 


a k c J_ n k n 
— O.OO ± U.oU 


u.uo 


a. oo _i_ n 71 
— o.oo ± U. / 1 


22 


9 ± 


.8 


Ujco oo 


onn/i a fK 
2UU4-Apr-Uo 


i omoAi/i 
1. /y /ZOj4 


t O/i _L n nc; 

— / .24 ± U.yo 


0.81 












2Uu4-May-U2 


28. /y5oUU3 


— o.oo ± 1.62 


1.22 


— 7.81 ± 1.63 


25 


9 ± 1 


2 


Uoco O / 


onn/i a n^ 
2UU4-Apr-Uo 


i "7 1 o oAon 


V (17 _L n Q.1 

— 7.U7 ± U.Of 


U.44 












2UU4-May-U2 


28. /l l J /yy 


— 5.70 ± 0.5s 


0.30 


/^ on 1 n oo 

— 6.38 ± 0.88 


18 


4 ± 


.4 


UOCO IJ 


onn/i Anr n/i 
2UU4-Apr-U4 


U.8o4U3 /O 


— O.oU ± U. /U 


n ?n 
U.oU 












2004-May-07 


33.6065432 


-6.16 ± 0.78 


0.60 


-6.48 ± 1.05 


55. 


6±3 


.0 


USco 100 


2004-Apr-06 


1.8179138 


-7.37 ±0.67 


0.44 












2004-May-02 


28.7752928 


-7.47 ±0.62 


0.36 


-7.42 ±0.91 


43. 


7±3 


.2 


USco 101 


2004-Apr-04 


0.8120734 


-3.25 ± 0.66 


0.44 












2004-May-02 


28.6591660 


-5.03 ± 0.60 


0.32 


-4.14 ±0.89 


19. 


1±0 


.3 


USco 104 


2004-Apr-04 


0.7850385 


-7.05 ± 0.61 


0.35 












2004-May-02 


28.6349480 


-7.62 ± 0.57 


0.27 


-7.33 ±0.83 


16. 


7±0 


1 


USco 109 


2004-Apr-05 


1.7453989 


-4.72 ± 0.50 


0.07 












2004-May-07 


33.6304878 


-4.70 ±0.51 


0.09 


-4.71 ±0.72 


8, 


6±1 


.2 


USco 112 


2004-Apr-04 


0.8552168 


-2.95 ± 0.50 


0.05 












2004-May-07 


33.5826368 


-3.41 ± 0.51 


0.07 


-3.18 ±0.71 


5. 


8±1 


2 


USco 121 


2004-Apr-24 


20.7059408 


-40.76 ± 0.90 


0.74 












2004-May-02 


28.6969641 


-42.48 ± 0.54 


0.20 


-41.62 ± 1.04 


17. 


6±1 


.3 


USco 128 


2004-May-13 


39.7978276 


-7.03 ± 0.51 


0.09 












2()04-May-17 


43.6108024 


-6.32 ± 0.51 


0.10 


-6.68 ±0.72 


3, 


6±1 


.1 


USco 130 


2004-May-09 


35.7724090 


-3.92 ± 0.94 


0.80 












2004-May-13 


39.8538830 


-3.90 ±0.67 


0.45 


-3.91 ± 1.16 


15. 


2±1 


1 


USco 132 


2004-May-13 


39.8268683 


-6.68 ± 0.53 


0.18 












2004-May-17 


43.6391138 


-6.61 ± 0.54 


0.20 


-6.64 ±0.76 


9. 


1±0 


.7 



100 times, and compute the standard deviation of the 100 RV val- 
ues. This should provide a good estimate of uncertainties in RVs 
cased by the uncertainties in the flux levels. We also estimated the 
RV uncertainties using the standard deviations on the mean of the 
two independent RV measurements obtaine d from the two co nsec- 
utive spectra taken on the same nig ht (c.f. Ijoergend Eiobh. We 
found that these estimates agreed well with those obtained from 
the Monte Carlo method. The 'average' radial velocities (RV) of 
the two epochs are also given in Tablef2| 

For each object and for each RV measurement, the deviations 
(ARV) from the average RV are computed and summarized in 
Fig.Qalong with their uncertainties (orrv) in order to aid the iden- 
tification of multiplicity candidates. Note that in computing ARV 
we do not require the knowledge of absolute/heliocentric radial ve- 
locities, but only the relative velocities (with respect to a template). 
We also cross-correlated object spectra from the two epochs (in- 
stead of using the spectra of the template LHS 49) in order to find 
ARV. The resulting cross-correlation functions were found to be 



too noisy (due to the relatively low S/N in the object spectra) for 
the peak positions of the cross-correlation curves to be reliably lo- 
cated. 

To identify an object with a RV variation with a statistical 
significance from our sa mple, we apply the method described by 
iMaxted & Jeffries! <2005l) . which we briefly summarize next, to our 
data. There are three steps in this method: (1) compute the \ 2 by 
fitting the two-epoch RV data for each object with a constant func- 
tion (a zero-th order polynomial), (2) compute the corresponding 
X 2 probability (p), (3) designate the object as a non-constant RV 
object or a binary candidate if p < 10 -3 (0.1 per cent). When 
computing x 2 , we use the uncertainties in the relative RV (ctrrv). 
The number of degree of freedom in the fitting procedure is obvi- 
ously 1. A si milar method was also u sed in a recent RV survey of 
VLM stars bv lBasri & Reinersl feOOfj) . 

We have computed p for all the objects (Table fJJ, and have 
plotted the results as a histogram of — logp, shown in Fig. [2] (ex- 
cluding the non-member USco 121; see explanation later). The fig- 
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Figure 1. Relative radial velocities (RVs) of objects measured in two different epochs. The vertical axes indicate the amount of deviation (ARV) from the 
'average' radial velocity (RV) in Table l2l and the horizontal axes indicate the time of the observation in heliocentric Julian date (HJD). The objects are 
considered to have a non-constant RV when the error bars of two data points do not overlap each other. The non-constant RV objects are considered as 
binary/multiple candidates. 



ure clearly shows two distinctive populations: one on the left (with 
small — log p values) occupied by the RV constant objects, and one 
on the right (with the large — log p values) occupied by the RV 
variable candidates. The expected distribution of p computed con- 
sistently with our uncertainty measurements is also shown in the 
same figure. The RV constant population on the left side reasonably 
matches the expected curve. To quantify this point, we compute the 
X 2 probability to test the goodness of the fit of the expected distri- 
bution to the histogram. We restrict the fit to the first 4 bins from the 
left (0 < — log 10 p < 2) of the histogram since the objects in the 
second bin from the right (2.5 < — log 10 p < 3) are 'potentially' 
RV-variable, although we have flagged them as the RV-constant ob- 
jects. In this analysis, we find the probability of 0.5 which indicates 
that the fit is reasonable. Further, if we re-normalized the expected 
distribution, accounting for that fact we neglected the objects in 
the second bin from the right, we find the \ 2 probability of the fit 
increases to 0.9. This, in turn, indicates that our uncertainties in 
RV values are reasonable, and an additional systematic error is not 
necessary in this analysis. This assertion is supported by the re- 



analysis of BD RVs, obtained by Joereens 1 2006b) using the same 
instruments setup as ours, bv lMaxted & Jeffrieslfeoolh . 

The boundary (p = 10 -3 ) between the RV variable and the 
RV constants se ems somewhat arbitrary, b ut here we simply adopt 
the definition of iMaxted & JeffrieU 120051) . Based on this criteria, 
there are four objects which show significant RV variations (out of 
17 samples) as one can see from the figure. They are USco 112, 
USco 128, USco 40, and USco 55, and considered as our prelimi- 
nary binary/multiple candidates. We will discuss the binary fraction 
and the expected binary detection probability later in Section|4| 

Finally, the histogram of RV for the objects in UScoOB is 
given in Fig. [3] The total number of the objects is 14. Note that 
USco 121 is excluded from the graph since it is identified as a 
non-member of the UScoOB associa tion based on the RV value 
(see Table |2j- Muzerolle et alj 1200 3T) also found it to be a likely 
non-member based on the radial velocity and the low lithium abun- 
dance. The distribution of the RVs in the figure was fitted by a 
gaussian function. We found that the standard deviation and the 
peak position of the radial velocity distribution are 1.0 km s -1 
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Figure 1. continued 



and — 6.3kms 1 respectively. The former is very similar to the 
standard deviation (0.9 km s _1 ) of the radial v elocity distributio n 
of 9 BDs and VLM objects in Cha I found bv ljoergensl <2006J) . 
They also studied the radial velocity distribution of more mas- 
sive 25 T Tauri stars in Cha I, and found the standard deviations 
(1.3 km s _1 ) is not significantly different from that of the brown 
dwarfs and the very low-mass objects. Unfortunately, we do not 
have the radial velocity measurements of the higher mass counter 
parts (T Tauri stars) in Upper Sco OB association. This is planned 
for future investigation since this is important for the study of the 
mass dependency of the kinematics in a young stellar cluster. 



According to the h ydrodynam i cal sim ulations of a low-mass 
star-forming cluster of iBate et alj 1200 3h which yields a stellar 
density of ~ 10 3 stars pc -3 , the rms dispersion (1-D) of the 
stars and the BDs is 1.2 km s~ . Similarly for the model with a 
higher stella r density (~ 10 4 star s pc~ 3 ), the rms dispersion is 
2.5 km s" 1 <Bate & Bonnell2005l) . The standard deviation of RV 
(1.0 km s" 1 ) found in our analysis is more comparable the lower 
stellar density model. 



3.2 Rotational velocities 

The rotational velocities of the objects were determined by measur- 
ing the widths of the cross correlation functions of the target spec- 
tra against a template spectrum from an object which is known to 
have a very small rotational velocity. The line broadening of the tar- 
gets is assumed to be dominated by rotational broadening. As in the 
cases for the radial velocity measurements, LHS 49 is chosen as the 
templ ate. Using its rotational period (P ~ 83 d. lBenedict & et alj 
1 19981) and radius (R* ~ 0.145 from the VLTI measurement 
bv lSeeransan et all2003l) . the rotational velocity of LHS 49 is esti- 
mated as v sini = 2irR*/P ~ 0.1 km s~ ; negligibly small. 

The width of the cross-correlation curves (cccf) are cal- 
ibrated with the rotational velocities (v sini) by cross corre- 
lating the template spectra against the same template spec- 
tra with ad ded rotation (convolv e d with a given v sinz ), as 
done by e.g.lTinnev & Reid ll998l) . lMohantv & BasrH <2003l) and 
I White & Basri 120031) . A linear limb-darkening law with a solar- 
like parameter (e = 0.6) was assum e d in th e formulation of the 
rotational profile described by iGravl 1 19921) . his Eq. 17.12. For 
each object, two measurements of rotational velocities are com- 
puted from two independent spectra obtained at different epochs. 
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Figure 2. Histogram of the x 2 probability (p) for fitting the observed (rel- 
ative) RV values with a constant (a horizontal line). The constant used in 
the fit is determined from the weighted mean of the two RV measurements 
for each object. The objects with p < 10 -3 are identified as non-constant 
or multiple, which appear in the right-most bin in the histogram. There are 
four objects which satisfy this condition. The expected distribution (dashed) 
for RV constant objects is also shown for a comparison, and is normalized 
to the total number of objects that have p > 10~ 3 . The match between 
the expected distribution and the histogram is reasonable (c.f. Section lTTT . 
indicating that our uncertainty estimates in RV values are also reasonable. 



Table 3. The x 2 probabilities (p) of each object being a RV constant, listed 
in ascending order of p. The ellipses represent the boundary (p = 10 — 3 ) 
between RV constant objects and RV variable objects. 
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As for the RV measurements, the mean and the standard deviation 
of the mean are used as the final rotational velocity and its uncer- 
tainty. The final results are recorded in Table|2| In general, our mea- 
su rements are in g o od ag reeme nt with the earlier measurements 
of iMuzerolle et all 1200 3l) and lMohantv 1 Jay_awardhana, & Basril 
120051) . given in Table Q] For example. IMuzerolle et alj 120031) 
found v sini = 16.8 ± 2.7 kms -1 for GY 5 while we found 
v sin i = 16.5 ± 0.6 kms -1 . 

The range of v sini found among our objects is 3.6- 



Figure 3. Histogram of the average heliocentric radial velocities of 14 US- 
coOB BD and VLM objects listed in Tablel2l(excluding USco 121, a non- 
member). The gaussian fit (dashed) of the radial velocity distribution gives 
a standard deviation of 1.0 km s -1 and the peak position of —6.3 km s -1 . 

55.6 k m s _1 , and a similar range is also found bv lMohantv et alj 
J2005l) . Fig. [4] shows the histogram of v sini distribution for the 
UScoOB objects (14 objects excluding USco 121, non member). 
The log-normal fit of this distribution gives the peak position at 
16.9 kms -1 with a standard deviation a = 27 .8 kms" 1 . Us- 
ing the v sini data in ljoergens & GuentheJ 1200 ll) . the same his- 
togram bin size used for UScoOB objects and a log-normal fit, we 
find the v sini distribution of the BD and VLM stars (8 objects) 
in Cha I peaks at 15.4 km s , and has the standard deviation of 
8.0 km s _ . The peak of the distribution is similar to that of Up- 
per Sco objects, but the standard deviation of the distribution is 
much smaller than that of our Upper Sco objects. The difference 
maybe due to the very small sample. A similar fit was applied to 
the v sin i dist ribution of 14 higher mass T Tauri stars in Cha I us- 
ing the data of Ijoergens & GuentheJ 1200 ll) . and we found a peak 
at 17.0 km s -1 with a standard deviation 25.9 km s _1 which are 
very similar to those of the Upper Sco brown dwarf candidates and 
VLM stars. 



4 BINARY FRACTION 

In Section ITT] we found 4 out of 17 (excluding USco 121; non- 
member) objects show a statistically significant RV variation, indi- 
cating that they are binary/multiple candidates. In order to estimate 
the uncertainty in the binary/multiple fraction fro m this relatively 
small sample, we will follow the method used bv lBasri & Reinera 
12006J) who considered the binomial distribution, Pb (x, n,p) of x 
positive event out of n trials with the probability p for a positive 
event in each trial. In our case, n = 17 (the number of sample) and 
x = 4. The peak of Pb (4, 17, p) curve suggests the binary fraction 
of 24 per cent, as its should be (4/17 = 0.24). The uncertainties 
was estimated by plotting Pb (4, 17, p), and finding the values of p 
at which Pb reduced to e -1 of the peak value. In this analysis, we 
find the binary/multiple fraction along with the uncertainties of our 
sample to be / = 24 jl^l per cent. 

Next, we investigate the range of semimajor axes of binaries 
(equivalently the range of binary periods) to which our RV sur- 
vey is sensitive. For this purpose, we will consider the detection 
probability for binaries or the RV variables, given the time sepa- 



8 R. Kurosawa et al. 



S 3 



0.1 - 



0.01 



10 20 30 40 50 60 

v sin i (km s ) 

Figure 4. Histogram of the rotational velocities (v sin i) of 14 UScoOB BD 
and VLM objects listed in Table l2l (excluding USco 121, a non-member). 
The log-normal fit (dashed) of the rotational velocity distribution gives a 
standard deviation 27.8 km s~ x and the peak position 16.9 km s — 1 . 

rations of two-epoch observations and the ranges of estimated pri- 
mary masses (c.f. Tables \T\ and The probability is calculated 
based on the simulated RV observations of binaries whose orbit are 
randomly selected from a mo del. A similar method w as used by 
iMaxted & Jeffries! J2005I) and iBasri & Reinersl fcfJOrJ) . The most 
important factors in determining the detection probability are the 
size of uncertainties in RV measurements (which we use the aver- 
age ctrrv from our observation), and the time separations of obser- 
vations. The smaller the errors in RV measurements, the larger the 
probability for a given binary orbit and a time separation of orbit. 
The larger the time separation of observations, the larger the upper 
limit of the semimajor axis to which an observation is sensitive. 
The average time separation of the two-epoch RV observations of 
our targets is 21.6 d. In the following, we will briefly discuss our 
model ass umptions and parameters which are essentially the same 
as those o flMaxted & Jeffrie! 120051) but with some simplifications. 

There are six basic parameters in our Monte Carlo simula- 
tion: primary mass (Mi), mass ratio (q), eccentricity (e), orbital 
phase (0), orbital inclination (i), longitude of periastron (A). The 
primary mass Mi is assigned from the adopted mass of the targets 
in TableQ with a uniform random deviation of ±0.002 Mq. The 
mass ratio is assumed to be uniformly distributed between q = 1.0 
and 0.2. The eccentricity e is a ssu med to be zero (circular orbits). 
Both lMaxted & Jeffries! <2005h and lBasri & Reinersl <200rj) found 
the detection probability is insensitive to the assumed distribution 
of q and e. The orbital phase is randomly chosen between and 1. 
The inclination i is randomly chosen from the cumulative distribu- 
tion of cosi. The longitude of periastron A is not necessary since 
we assumed e = 0. 

In order to compute the detection probability as a function of 
semimajor axis a, we take the following procedure: (1) for each 
object in our targets, we randomly select 10 5 binaries using the as- 
sumption stated above for a given value of a, and compute the RV 
of the primary (Vi), (2) evolve the orbit by the time separation of 
the RV measurements used in the observations for this object, and 
take another simulated measurement of RV (V2), (3) from Vi, V2 
and the average uncertainty in RV (ctrrv) from the observations, 
we compute the x 2 probability p, and flag the trial as a detection 
if p < 10 -3 as done for the real data, (4) find the fraction of de- 
tections out of all random trials, (5) repeat \-4 for the range of a 
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Figure 5. The binary detection probability (solid) as a function of the semi- 
major axis (a), based on Monte Carlo simulations of RV measurements, is 
shown. For each object in our target list, the simulation was performed us- 
ing the time separation (At) actually used in our two-epoch observations 
(c.f. Tablel2l and the estimated primary mass (c.f. TablelTl. The final detec- 
tion probability is obtained by averaging over the simulated observations of 
all the objects. The average time separation in the two-epoch observations 
is ~ 21.6 d, and the average uncertainty in RV is 0.42 km s —1 . The prob- 
ability remains fairly constant up to a ss 0.1 au, and it rapidly decreases 
beyond a ks 0.3 au. The 80 per cent detection probability is achieved for 
binaries with a < 0.1 au. 



between 10 -3 and 10 au, and (6) repeat 5 for the all targets, and 
find the detection probability averaged over all targets as a function 
of a. 

The result of the simulation is shown in Fig. B] The detection 
probab ility curve shown here is very similar that o flBasri & Reinersl 
(2006) for a constant time separation (20 d) case (see their Fig. 4). 
The probability remains fairly constant up to a ~ 0.1 au, and it 
rapidly decreases beyond a ~ 0.3 au. This turning point will in- 
crease if we had used larger time separations in our observation. 
From this figure, we find that the 80 per cent detection probability 
up to for binaries with a < 0.1 au, given our time separations and 
the average uncertainty (~ 0.42 km s ) in the RV measurements. 

An alternative explanation for the RV variations found in the 
binary candidates abo ve is a type of stellar pulsations found by 
IPalla & Baraffel l2005h . They studied the non-adiabatic, linear in- 
stability of very VLM stars and BD during the deuterium burning 
phase in the core, and found unstable fundamental modes in the 
time-scale between ~ 1 and ~ 5 h for the object mass from 0.02 
and 0.1 Mq. By using their pulsation periods and by assuming that 
the pulsation amplitudes is 5 per cent of the radius of the objects, 
we find that the expected RV variation can be > 6kms~ . As 
mentioned in Section|2| for each object at a given night, two sep- 
arate spectra are obtained consecutively (with time separations of 
~ 0.5 h.). We have measured the RV from each spectrum to check 
if there is a large jump in the RV measurements on a short time- 
scale. The size of the shifts in the RV values measured in the two 
consecutive spectra are found in between 0.02 ± 0.10 (USco 128) 
and 2.04 ± 1.67 km s _1 (USco 66), and the average size of the 
shifts 0.36 km s _1 , which is m uch smaller than the ex pected RV 
variation based on the models o f|Palla& Baraffel <2005l) : hence, we 
found little indication of the pulsations in our sample. 
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5 CONCLUSIONS 

We have presented two-epoch RV survey of 18 young BDs and 
VLM stars (0.02 M & < M* < 0.1 Mq) in UScoOB and p Oph 
dark cloud core using the high resolution UVES echelle spec- 
troscopy at VLT. The average time separation of RV measurements 
are 21.6 d, and our RV survey is sensitive to binaries with sepa- 
ration smaller than 0.1 au. One of our targets, USco 121, is most 
likely a non-member of the UScoOB association based on the devi- 
ation of the RV from the result of the po pulation in the assoc iation. 
A similar conclusion was found by iMuzerolle et alj l2003h from 
their RV study and the low lithium abundance. 

We found 4 (USco 112, USco 128, USco 40 and USco 55) out 
of 17 objects as our binary/multiple candidates. This corresponds 
to the binary fraction of 24^^ per cent for the binary separation 
a < 0.1 au. The recent high-resolution imaging survey of brown 
d warfs and very low- mass objects (M5.5-M7.5 ) in the UScoOB 
bv lKraus et alj|2005l) (which was not known to authors at the time 
of our observation: April-May, 2004) confirms that USco 55 and 
USco 66 are multiple systems, and USco 109 is most likely a mul- 
tiple system. Their projected separations are a > 4.0 au which is 
well above our sensitivity limit of 0. 1 au; therefore, it is not surpris- 
ing that we did not identify USco 66 and USco 109 as multiples. 
Interestingly, we also found USco 55 as a candidate for a multiple 
system. This may indicate that it is a triple system with one of the 
objects located within 0. 1 au of the primary. In addition, we iden- 
tify USco 112 and 128 as binary candidates, but they did not find 
them as multiples. Further they found USco 67, 75, 130 and 132 as 
non-multiple, as did we. 

We found the RV dispersion (1.0 km s" 1 ) of the objects in 
UScoOB is very similar to that of the BDs and VLM st ars in 
Chamaeleon I (Cha I) previous study by IjoergenJ <2006iJ) . The 
rotational velocities (v sini) of the samples were also measured. 
The distribution of v sini for the UScoOB objects peaks around 
16.9 km s -1 which is also similar to that of the Cha I population 
found bv ljoergensl l2006af) : however, the dispersion of v sini for 
the UScoOB objects (27.8 km s^ 1 ) is found to be much larger than 
that of the Cha I objects (8.0km s _1 ) 

Follow-up spectroscopic observations of the binary candi- 
dates presented here are planned in near future. There are only 
a few RV variable binary candidates identified in earlier survey s 
iGuenther & WuchterfcOoSlKenvon et all2005tljoergenj2006hl) . 
AT lBTIrga^ser^HuTl2006Tr points out most of the current RV and 
imaging surveys use samples from magnitude-limited survey, but 
one should attempt to use the samples from volume-limited sur- 
vey in order to find a correct statistics on binary parameters more 
straightforwardly, i.e. without correcting for bias. It is possible 
that our binary fraction is an overestimate since we are biased to 
brighter objects. However the lack of a double-line system sug- 
gests that a primary-to-secondary luminosity ratio of 3> 1 for the 
RV variables found in this paper. 
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